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Atomic force microscopy was used to investigate ferroelectric switching and nanoscale domain 
dynamics in epitaxial Pb(Zro.2Tio.8)03 thin films. Measurements of the writing time dependence 
of domain size reveal a two-step process in which nucleation is followed by radial domain growth. 
During this growth, the domain wall velocity exhibits a v oc exp —(l/E)^ dependence on the electric 
field, characteristic of a creep process. The domain wall motion was analyzed both in the context 
of stochastic nucleation in a periodic potential as well as the canonical creep motion of an elastic 
manifold in a disorder potential. The dimensionality of the films suggests that disorder is at the 
origin of the observed domain wall creep. To investigate the effects of changing the disorder in the 
films, defects were introduced during crystal growth (a-axis inclusions) or by heavy ion irradiation, 
producing films with planar and columnar defects, respectively. The presence of these defects was 
found to significantly decrease the creep exponent n, from 0.62 - 0.69 to 0.38 - 0.5 in the irradiated 
films and 0.19 - 0.31 in the films containing a-axis inclusions. 

PACS numbers: 77.80.Fm, 77.80.Dj, 68.37.Ps, 61.80.Jh 



I. INTRODUCTION 

Ferroelectric materials are widely used in modern tech- 
nologies in order to exploit their ferroelectric, piezoelec- 
tric and pyroelectric properties. Standard devices gener- 
ally use microscopic patterned electrodes on ferroelectric 
ceramics, single crystals or films to produce actuators, 
filters, resonators, sensors and memories. However, the 
continuing demand for miniaturization in acoustic ap- 
plications in ultra- high density information storage 
H,0J, and in micro-electro-mechanical systems (MEMS), 
has made alternative solutions necessary. In terms of 
material growth, techniques like rf-magnetron sputter- 
ing, pulsed laser deposition, and molecular beam epitaxy 
with in-situ characterization £| have been developed to 
precisely control the structure and thickness of materi- 
als. Today, these techniques allow epitaxial thin films 
with atomically flat surfaces and interfaces to be grown, 
both in all-oxide structures and, more recently, on top 
of silicon p|. Combined with these techniques, increased 
understanding and control of oxide growth modes has 
permitted self-assembly processes and nanostructuring 
in single crystals and thin films [jj, a novel approach 
to miniaturization. Low-temperature metallo-organic 
chemical vapor deposition, more easily adaptable to in- 
dustrial applications, has also been actively researched to 
produce uniform films over large areas with low fatigue 
0. As device dimensions become reduced in extremely 
miniaturized systems, an important issue is the ability to 
locally access and control ferroelectricity. From its incep- 
tion, there has been significant interest in using atomic 
force microscopy (AFM) [HUE! for this purpose, since 
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it allows control over ferroelectric domain structure at 
the nanoscopic scale required by ever smaller systems. 
For such systems, ferroelectric domain stability, domain 
dynamics in the presence of electric fields or temperature 
variations, and issues of domain growth are the main con- 
cerns. Therefore, understanding the fundamental physics 
of ferroelectric domains in thin films is of crucial impor- 
tance. 

In this respect, it is useful to consider domain walls 
in ferroelectric materials within the broader framework 
of elastic disordered systems in the presence of an exter- 
nal force. Indeed, ferroelectric materials are character- 
ized by energetically equivalent, degenerate ground states 
(tetragonal, in the case of the perovskite materials used 
in this study), separated by an energy barrier. In a fer- 
roelectric ground state the center of gravity of positive 
charge in the unit cell is displaced relative to the center 
of gravity of the negative charge, leading to a stable re- 
manent polarization reversible under an electric field, as 
shown schematically in Fig.Q] Regions of opposite polar- 
ization are separated by thin interfaces, or domain walls. 
The application of an electric field asymmetrizes the fer- 
roelectric double well potential, favoring one polarization 
state over the other by reducing the energy necessary 
to create a nucleus with polarization anti-parallel to the 
field, thus promoting domain wall motion. These domain 
walls can be considered as elastic objects, whose surface 
tends to be minimized in order to decrease the total wall 
energy. However, such elastic objects may meander from 
an optimal flat configuration in order to take advantage 
of particularly favorable regions of the potential land- 
scape, which can vary due to disorder in the ferroelectric 
film, as well as the commensurate potential of the crys- 
tal lattice itself. The behavior of such systems is thus 
governed by the competition between elasticity and the 
effects of pinning. Although at T = the domain wall is 
pinned by the disorder until a critical force f c is reached, 
at finite temperatures it can be driven by forces below 
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FIG. 1: Schematic diagram of a Pb(Zro.2Tio.g03 unit cell in 
two degenerate ground states separated by an energy barrier 
Uo- The Pb and Ti/Zr ions are positively charged and the 
O ions are negatively charged, giving opposite polarization 
directions Pdown and Pup. 



f c , since barriers to motion, however high, can always 
be passed via thermal activation. How the domain wall 
moves when driven at small forces is of theoretical and 
practical interest. A simple answer jllj is motion gov- 
erned by a linear response of the form v oc e _/3A / due to 
thermal activation above pinning barriers A, where (3 is 
the inverse temperature 1 /kpT and / <C / c is the driving 
force. However, subsequent detailed analysis showed that 
both periodic |12l Il3| and disordered pinning potentials 
(IH ITM IT3 | can lead to diverging barriers and thus to a 
nonlinear "creep" response, in which the velocity is, for 
small temperatures and for forces much smaller than the 
critical force, of the form v oc exp {—(3U c {f c /jY). U c is 
a barrier height and \i is an exponent characteristic of 
the mechanism responsible for the creep. In a periodic 
potential creep occurs only for interfaces of dimension 
d = 2 (sheets), and /i = 1 [13|. For creep in a random 
potential /u = (d — 2 + 2£) / (2 — £) , where d is the dimen- 
sionality of the interface and £ its equilibrium roughness 
exponent (see Section^ . Note that the nature of the dis- 
order only enters through the static roughness exponent 
£. Measurements of domain wall motion and equilibrium 
roughness configuration in ultrathin magnetic films have 
verified the creep law prediction for interfaces, with the 
measured exponent fi = 0.25 in good agreement with 
the expected theoretical values for this system ^tJ- In 
periodic vortex systems, precise determination of the ex- 
ponents has been complicated by the multiple scales of 
the problem 0, , although results in agreement with 
the theoretical predictions of ( = have been observed 
[T^|. High quality epitaxial ferroelectric thin films are 
therefore another interesting model system to study such 
phenomena, because their thickness and crystalline qual- 
ity can be precisely controlled, and because AFM allows 
nanometer resolution of their domain configuration. 



In this paper, we present a review of our studies of 
non-invasive, local AFM control of ferroelectric polar- 
ization in epitaxial Pb(Zro.2Tio.s)03 (PZT) thin films 
and its possible applications. We will focus specifically 
on the behavior of domain walls and its relation to do- 
main growth and stability. We expand on the discussion 
of our previous results demonstrating that domain wall 
motion in thin films is a disorder-controlled creep pro- 
cess. In particular, we discuss the differences between the 
usual stochastic nucleation scenario proposed by Miller 
and Weinreich to explain domain wall motion observed 
in bulk ferroelectrics [l^, and the creep scenario due 
to the competition between a disorder potential and the 
elasticity of the wall. Finally, we report on new studies 
of ferroelectric materials in which domain dynamics are 
measured in the presence of macroscopic defects, both 
columnar, in the form of heavy ion irradiation tracks, 
and planar, in the form of a-axis inclusions in the c-axis 
oriented films. We find that these defects, macroscopic 
in relation to the thickness of the domain wall, nonethe- 
less have a significant effect on domain wall dynamics, 
lowering the values of the critical exponent for domain 
wall creep from 0.62 - 0.69 to 0.38 - 0.5 in irradiated 
films, and down to 0.19 - 0.31 in films containing a-axis 
inclusions ppf. 



II. THIN FILM FABRICATION AND AFM 
WRITING 

The ferroelectric materials used in this study were c- 
axis oriented PZT thin films, epitaxially grown by off- 
axis rf-magnetron sputtering onto conductive (0.5 % wt) 
Nb:SrTi0 3 (001) substrates in an Ar:0 2 flow at 180 
mTorr, and at substrate temperatures of ~500°C. Mul- 
tiple samples with thicknesses varying from 29.0 to 130.0 
nm were grown. X-ray characterization of the films, such 
as the 9-29 diffractogram shown in Fig. |2 revealed high 
crystalline quality, with (f> scans (not shown) confirming 
epitaxial "cube-on-cube" growth of the ferroelectric ma- 
terial on the substrate. The multiple orders of satellite 
reflections around the principal 001 PZT peak are due 
to the finite size of the sample, and allow us to precisely 
determine the thickness of the films. Measurements of 
sample topography showed uniform surfaces with root- 
mean-square (rms) roughness of 0.2 - 0.3 nm over 5x5 
fxm areas as shown in Fig. [31 where the vertical scale is 
equal to 1 % of the 54 nm thickness of the film itself. High 
crystalline quality and smooth surfaces are advantageous 
for effective writing and imaging of ferroelectric domains, 
since the presence of morphological defects could perturb 
the tip-sample interaction and thereby complicate the use 
of the films as model systems for domain wall studies. 

To locally control the ferroelectric polarization state, 
a metallic AFM tip was used as a mobile top electrode, 
while the conductive substrate functioned as a bottom 
electrode. Applying a voltage between the tip and the 
substrate results in a local electric field across the fer- 
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FIG. 2: 8-26 scan of the PZT 001 peak in a 50 nm epitaxial 
thin film. Multiple orders of finite size reflections allow precise 
measurement of film thickness. 




FIG. 3: AFM topography of a 54 nm PZT film showing 0.22 
nm rms roughness over a 5 x 5 /im 2 area. 



roelectric thin film, favoring the polarization direction 
parallel to the field. This leads to polarization switching 
if the field is high enough (or applied for long times), as 
shown schematically in Fig. 0] To image the resulting do- 
mains the phase contrast of the local piezoresponse of the 
film, excited by the application of a small oscillating volt- 
age via the AFM tip, (detailed in Fig.JSJ is measured by 
piezo- force microscopy (PFM). All measurements were 
performed on commercially available Veeco Multimode 
equipment. 

Large areas, essentially limited only by the scan area 
of the AFM and ultimately, the size of the sample it- 
self, can be polarized by applying a constant voltage to 
a tip scanning in contact with the surface 9]. One can 
also alternate positive and negative voltages applied for 
a fixed duration with respect to the scanning speed of 
the sample to create arrays of lines with opposite polar- 
ization, such as those shown in Fig. HJ1 where the line 
width is 925 ± 15 nm. Such linear domain structures 
have been used to develop a new type of interdigitated 
transducer for a pro totype high-frequency surface acous- 
tic wave device [21j. Frequencies up to 3.4 GHz have 




FIG. 4: Schematic representation of local polarization switch- 
ing in a ferroelectric thin film. The application of voltage 
pulses of magnitude AV and duration At between a metallic 
AFM tip and the conducting substrate produces nanoscale 
circular domains penetrating through the thickness of the 
film. 




FIG. 5: Schematic representation of the two antisymmetric 
polarization states in the tetragonal ferroelectric PZT unit 
cell. Applying a small AC voltage with an AFM tip excites a 
local piezoresponse, detected using a lock-in technique. The 
two polarization states (Pup and Pdown-) respond 180° out 
of phase with each other, thus allowing a phase contrast image 
of the ferroelectric domains to be obtained. 

been demonstrated, but higher frequencies are in prin- 
ciple readily accessible since AFM writing allows very 
small line widths. For instance, using carbon nanotube 
bundles attached to a metallic AFM tip to achieve a very 
high aspect ratio, we were able to write a network of 
lines as small as 20-30 nm |22^. Small circular domains 
can also be created by applying short voltage pulses to 
a stationary tip in contact with the sample, at desired 
positions in a uniformly prepolarized area. As we have 
previously shown, such domains can be written in arrays 
with densities as high as ~30Gbit/cm 2 , with each domain 
(bit) in such an array individually accessible and fully 
reversible under subsequent voltage pulses [2^, an im- 
portant consideration from the point of view of dynamic 
non-volatile memory applications. Subsequently, Cho et 
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FIG. 6: PFM image of AFM- written linear ferroelectric do- 
mains 925 ± 15 nm wide. Similar domain structures have 
been used in a prototype surface acoustic wave device with 
GHz frequencies. 



al. have demonstrated individual domains ~15 nm in 
diameter in ferroelectric single crystals which would give 
densities as high as 0.5 Tbit/cm 2 [24J when projected into 
a standard array. By identifying the parameters control- 
ling domain size as the writing time and writing voltage 
(the duration and magnitude, respectively, of the voltage 
pulse applied to create the domain) , and the confinement 
of the electric field at the AFM tip our studies address an 
active area of interest for ferroelectric non- volatile mem- 
ories: how to achieve the required ultrahigh information 
densities. 

Two other key issues for information storage applica- 
tions are the minimum switching time and the stability 
of the resulting domains. Using the AFM, we have writ- 
ten domains with pulses as short as 5 ns, with radii of 
18.5 ± 4 nm. This is not a fundamental limit, however, 
but rather an experimental consideration related to the 
minimum pulse time achieved by our pulse generator and 
the RC characteristics of our system, which may lead to 
distortion of the applied pulse at very short times. Recent 
studies by Li et al. using a photoconductive switch and 
femtosecond laser illumination suggest that the actual 
switching time for ferroelectric domains may be as short 
as ^70 - 90 ps [25]. Our studies have also consistently 
shown very high stability of AFM-written ferroelectric 
domains. An array of 16 domains written with 500 ns 
pulses was imaged immediately after writing (Fig.[7fa)), 
showing well defined homogeneous domains with regular 
spacing and radii of 29.2 ±4 nm. The array was followed 
with measurements performed at 1 week (Fig. 0b)), 2 
weeks (Fig. [7fc)) and 1 month (Fig. [7^d)), with radii of 
27.0 ±4 nm, 27.2 ±4 nm and 27.7 ±4 nm respectively. 
No change, backswitching, or spontaneous disappearance 
of the domains were observed. Additionally, the arrays 
of lines used for surface acoustic wave devices were mea- 
sured over a period of 4 months, after lithographic pat- 
terning, etching of areas outside the filter region, and the 
application of voltage signals larger than those used for 
the piezoelectric response in order to excite surface acous- 
tic waves. Once again, the domain structures remained 
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FIG. 7: PFM images of the same array written with 500 ns, 
12 V pulses, straight after writing (a), 1 week later (b), after 
2 weeks(c) and at the end of the experiment, a month later 
(d). The average domain radii are 29.2 nm, 27.0 nm, 27.2 nm 
and 27.7 nm respectively, unchanged within the ±4 nm error 
of the experiment. 

completely stable [2l]. Moreover, similar linear domains 
heated to 440° C did not show any change when imaged 
after subsequent cooling j2(J. It is therefore likely that 
pre-written ferroelectric domain structures in thin films 
can subsequently be incorporated into further processing 
steps up to, and possibly well beyond this temperature. 
As shall be explained in more detail, such high stabil- 
ity in zero applied field is in fact inherent to the glassy 
behavior of an elastic disordered system. 

III. EXPERIMENTAL OBSERVATION OF 
DOMAIN WALL CREEP 

The local control of ferrolectric polarization provided 
by AFM, together with its nanometer resolution, also 
make it a powerful tool for fundamental studies of do- 
main dynamics and dependence on writing parameters 
such as the writing time. As previously reported |23| . we 
observed a strong dependence of domain radius on the 
writing time for times longer than ~20 fis. For shorter 
times, domain radii appeared to saturate at ~20 nm, a 
size we relate to the 25 - 50 nm nominal radius of cur- 
vature of the AFM tip used for these experiments. For 
each writing time, 16- or 25-domain arrays were writ- 
ten with 12 V pulses, and the domain radius was calcu- 
lated by averaging over the measured domain radii along 
their vertical and horizontal axes, with a rms error of 
~10%. All imaged domains appeared homogeneous and 
well defined and no randomly nucleated domains were 
observed. These data suggest a two-step switching pro- 
cess: first, rapid nucleation and forward growth across 
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the thickness of the sample occur under the AFM tip; this 
event is followed by slower radial motion of the domain 
wall outwards, perpendicular to the polarization direc- 
tion, in order to increase domain size. We analyzed this 
radial domain wall motion in the framework of a pinned 
elastic system by comparing the velocity and the driving 
force exerted on the wall, in our case due to the electric 
field E applied by the tip. We considered arrays written 
with consecutive pulse durations, extracting the domain 
wall velocity as v = r ^^~^ tl - > , the difference in domain 
radii at the two subsequent writing times divided by the 
difference in the writing times themselves. The electric 
field distribution was obtained by modelling the tip as 
a charged sphere, with radius a taken as equal to the 
domain saturation size of 20 nm. Applying a voltage V 
to the tip at the surface of the ferroelectric film with di- 
electric constant e produces a charge q = AireeoaV on the 
model tip. Taking into account both the effect of the film 
and the conductive substrate, we are able to find the field 
Ei_{r, z) at any point (r, z) within the film, r is the hor- 
izontal distance (in the plane of the film) away from the 
center of the spherical tip, and z is the depth within the 
film (up to thickness A) from the center of the tip. In our 
experiments, the domain radii remain comparable to the 
size of the tip, so further simplification can be obtained 
by considering only the first order of image charge reflec- 
tions in the film and the substrate. Since the voltage drop 
V across the film is simply the integral of this field over 
the film thickness V — Jq +X E±(r, z)dz we can define the 
average field across the film E(r) — E±(r), which shows a 
1 jr dependence in this first order approximation. As one 
moves further away from the tip, a crossover to higher or- 
ders of r dependence in the denominator is expected for 
the field. This simplified model shows reasonable agree- 
ment with a more accurate numerical simulation of a hy- 
perbolic tip in contact with a ferroelectric film grown on 
a metallic substrate. Here, a 1/r dependence of the elec- 
tric field is observed out to r ~ 150 nm, for a tip 1 /im 
high with a = 20 nm radius of curvature. Beyond this, 
a cross over to 1/r 2 and even steeper field decay occurs. 

Although the field E(r) is highly inhomogencous at 
large length scales, it can be taken as constant over the 
very small thickness (on the order of a lattice spacing) of 
the ferroelectric domain wall. One can thus relate the ve- 
locity v(r) of the domain wall at a distance r from the tip 
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to a field E(r) = where 
nius plot of the velocity against the inverse field, shown 
in Fig. 03 reveals that our data are in good agreement 
with a creep behavior 
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over multiple decades of velocity, from 10 -3 to 10~ 9 m/s, 
and for fields varying from 10 7 to 5 x 10 8 V/m 48]. The 
values of the dynamical exponent \x were found to be close 
to 1 for the original three samples (2?|], and closer to 0.7 
for seven later samples - grown and measured under the 
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FIG. 8: Domain wall speed as a function of the inverse applied 
electric field for 37.0, 54.5, 81.0 and 95.4 nm thick films. The 
data agree with the creep equation v ~ exp [— j^rp (■^■) M ] 
with /i = 0.93, 0.62, 1.19, and 0.70 respectively. Fits of the 
data to logw = A(1/S) M are shown for the 54.5 and 95.4 nm 
films. 



same conditions. These data were the first indication 
that domain wall motion in ferroelectric thin films was 
a creep process, and led us to investigate its microscopic 
origins. 



IV. DOMAIN WALL CREEP IN A 
COMMENSURATE POTENTIAL 



Early studies [28J, [29j of domain growth carried out by 
optical and etching techniques on bulk samples reported 
a non-linear electric field dependence of the velocity v cx 
exp— 1/E known as Merz's law (with implicit values of 
yu. = 1, if these results are to be considered in the general 
framework of creep). At the time, a phenomenological 
theory based on the stochastic nucleation of new domains 
at existing domain boundaries was put forward by Miller 
and Weinreich to explain the observed behavior [12j . The 
wall moves forward due to the formation of a nucleus as 
shown in Fig. EI The energy change due to the formation 
of a nucleus is 



AF — — 2P S EV + a w A + [/depolarization 



(2) 



Nucleation would occur and the domain wall would move 
if the energy gain, due to switching a volume V of fer- 
roelectric with spontaneous polarization P s to the polar- 
ization state energetically favorable with respect to the 
direction of the applied field E, would balance the en- 
ergy cost of extending the surface A of the domain wall, 
with a surface energy density of a w , as well as the in- 
curred depolarization energy cost V depolarization- In fact, 
this mechanism is identical to the one of an elastic 
manifold weakly driven in a periodic pinning potential 
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FIG. 9: Schematic drawing of a triangular step domain on a 
180° domain wall, as described by Miller and Weinreich ^3 
The applied electric field E is parallel to Pdown. 



(tilted washboard potential, as described for example in 
[l3j). The nucleus thickness (c in Fig. I^J is the distance 
between two mimima of the periodic potential given by 
the lattice spacing of the ferroelectric crystal. 

For small electric fields (E — ► 0) a large nucleus can be 
expected since the energy gain due to the displacement 
of a nucleus into the neighboring pinning valley grows 
with the volume of the nucleus, while the energy cost 
essentially scales with its surface. Since V ~ L d while 
A ~ L 1 , where L is the extension of the nucleus, and d 
the dimensionality of the elastic interface, two cases oc- 
cur. For a one-dimensional manifold (string, d = 1), the 
nucleus consists of two point-like kinks, whose activation 
energy therefore always remains finite, and the system 
exhibits a linear response under small driving forces. For 
a two dimensional manifold, on the other hand, mini- 
mizing gives L* r~j 1/E, showing that the size of the 
nucleus grows as the electric field decreases. The en- 
ergy barriers the nucleaus has to overcome thus grow 
as A(E) ~ 1/E, using (J2J, giving a non-linear response 
with v cx exp—l/E. The stochastic nucleation scenario 
proposed by Miller and Weinreich can thus explain the 
observed non-linear response of the domain wall only if 
the domain wall itself is a two-dimensional surface em- 
bedded in a three-dimensional crystal. This means that 
the dimensions of the nucleus, at a given field E, have to 
be smaller than the thickness of the system. Otherwise, 
the energy of the nucleus saturates and the one dimen- 
sional wall case and a linear response are recovered. It 
is also important to note that if the creep consists of 
motion in a periodic potential the creep exponent is con- 
strained to be n = 1. As already mentioned, this par- 
ticular scenario is microscopically related to the intrinsic 
periodic potential of the ferroelectric crystal itself acting 
to pin the domain wall. The strength of this potential 
was calculated in ab-initio studies of 180° domain walls 
in PbTiOa, showing that the wall energy varies from 132 
mJ/m 2 to 169 mJ/m 2 depending on whether the domain 
wall is centered on a Pb-0 or T1-O2 plane in the crystal 



[3fJ. The influence of such a periodic potential is due to 
the extreme thinness of the domain wall in ferroelectric 
materials, in contrast to magnetic systems where the do- 
main wall is much larger than the atomic length scales 

m 

In order to test whether the observed creep behav- 
ior is indeed due to the nucleation process, we calcu- 
lated the size of the critical nucleus, following the for- 
mulation derived by Miller and Weinreich for the ener- 
getically most favorable dagger-shaped nucleus of hori- 
zontal extension a, height I and thickness c forming at 
an existing 180° domain wall, as shown in Fig EI 12], 
where P s is the polarization, b the in plane lattice con- 
stant, e the dielectric constant of PZT at ambient con- 
ditions, and E the applied electric field. The depolar- 

ization energy can be written as {/depolarization 



with cr p = (4P 2 61n(0.7358a/6))/e 



By minimizing 



the free energy change due to nucleation with respect 
to the dimensions of the nucleus a and I, with c taken as 
equivalent to the lattice constant b (the distance between 
two minima in the periodic crystalline potential) , the size 
of the critical nucleus a* and I*, as well as the activation 
energy AF*, can be calculated as |T^ |: 
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To compute the actual values, we take the standard pa- 
rameters for PZT (P s = 0.40 C/m 2 , e = 100, b = 3.96 
A ), and the ab-initio value for the domain wall energy 
density 50] a w — 0.132 mJ/m 2 . In our case, the applied 
electric field varied from ~ 2 to 20 MV/m (with the factor 
10 correction), depending on the thickness of the sample 
used and the distance from the AFM tip, with the most 
intense fields for thin films and small domains. Corre- 
sponding values of o~ p were between 1.6 and 0.9 J/m 2 . 
Since a p is therefore greater than a w , following Miller 
and Weinreich, the expressions for the critical values can 
be simplified to: 



AF* 



3 P S E 
„ 1/2 1/2 

V3P S E 
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For the field range used, these equations would give crit- 
ical values of a* ~ 12.5 - 125 nm and I* ~ 53 - 710 nm. 
These results imply that for the given value of field, the 
vertical size of the critical nucleus would exceed the thick- 
ness of the film itself. This suggests that the films are in 
a 2-dimensional limit, with the domain walls acting as a 
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FIG. 10: Domain wall as elastic manifold trapped in a random 
disorder potential. In equilibrium, the domain wall exhibits a 
characte ristic roughness, mea sured by the correlation function 
B(z) — {[u(z + L) — u(z)] 2 ) of the displacements u(z) from 
elastically ideal flat configuration with respect to the length 
L of domain wall. In the presence of a small driving force 
/ < f c due to the applied electric field E, the domain wall 
will move to the next favorable configuration in the disorder 
potential, as shown by the dotted line, via a glassy creep 
motion. 



quasi-one-dimensional manifold, for which the Miller and 
Weinreich stochastic nucleation model, or alternatively, 
weakly driven motion through a periodic potential, could 
not explain the non-linear response observed. Finally, the 
values of the dynamical exponent we observe, generally 
not equal to one, are also a strong indication that an al- 
ternative microscopic mechanism for the observed creep 
process should be considered. 



DOMAIN WALL CREEP IN A RANDOM 
POTENTIAL 



In the alternative scenario of a canonical "glassy" sys- 
tem, an elastic manifold is weakly collectively pinned by 
the quenched disorder potential present in the medium, 
with important consequences for both its static and dy- 
namic behavior. Since disorder is always present in any 
realistic system it is a feasible source of the observed 
domain wall creep in ferroelectric thin films. Both for 
vortices in superconductors 01 an d in magnetic systems 
|17| , disorder is clearly the driving mechanism behind the 
observed creep behavior. In PZT, a solid solution of 20% 
PbZrC>3 in 80% PbTiC>3, the presence of Zr atoms is a 
possible source of disorder, although preliminary studies 
of domain wall dynamics in pure PbTiC>3 show similar 
behaviour to that observed in PZT. Vacancies and other 
defects in the lattice structure are also likely sources of 
disorder. In ferroelectric films the presence of disorder 
would dominate domain wall behaviour for both one and 
2 dimensional walls at large scales. However, given the 
thinness of the domain wall, we note that the periodic 
potential of the crystal is also present in the problem. In 
principle, ferroelectric films could therefore ultimately be 
used to study the competition between a periodic poten- 
tial and disorder (see e.g. |3lj')- 



In order to analyze the effects of disorder on domain 
wall motion let us again consider the energy for a segment 
of ferroelectric domain wall of length L displaced by u(z) 
from the elastically ideal flat configuration (where z are 
the internal coordinates, such that the total spatial di- 
mension d = z + 1), as shown on Fig. 1101 The energy 
scales as HlJ 



U(U,L) = a w U 2 L d - 2 - (JdisordcrM - 2P s EL a 



(5) 



where the first term describes the elastic e nerg y contri- 
bution, and is expressed for a local elasticity |§2 • A more 
accurate description of long range forces, such as dipolar 
forces, modifies the elasticity and amounts to replacing d 
by (3d — l)/2 in the following formulas(see e.g. [3l| and 
ref. therein). The second term is due to pinning by the 
disorder potential, and the third is the energy due to the 
application of an external electric field, f/dis. depends on 
the precise nature of the disorder. For example for "ran- 
dom bond" disorder, equivalent to defects which locally 
modify the ferroelectric double well potential depth, one 
can model the disorder by a random potential acting at 
the position of the interface and 



C^dis 



rdcr 



d d zV(u(z),z) 



(6) 



Another form of disorder, the so called "random field", 
occurs when defects locally asymmetrize the ferroelec- 
tric double well potential, leading to a different form for 
-Edis.- If the disorder is weak, the central limit theorem al- 
lows its approximation by a Gaussian random potential. 
The disorder is then only characterized by its correlation 
length rt and the strength of the random potential. In 
the absence of an external electric field, the configuration 
of the domain wall results from the competition between 
the elastic forces and the random potential. This configu- 
ration can be characterized by measuring the correlation 
function of relative displacements 



B{L) = ([u(z + L) u(z)f) = e 



2C 



(7) 



where (• • • ) denotes thermal averaging (thermodynamic 
equilibrium) and 7TT denotes an ensemble average over 
the realization of the disorder. In a realistic experimen- 
tal situation the ensemble average is performed by aver- 
aging over all pairs of points separated by a distance L, 
assuming that the system is self-averaging. B(L) shows a 
power-law growth with different exponents. For r smaller 
than a characteristic length, the Larkin length L c |32ll33| . 
B(L) grows as B(L) ~ L A ~ d . Below this length there is 
no metastability and no pinning of the elastic interface. 
Above the Larkin length, the growth still follows a power 
law, but with an exponent 2£ (B(L) ~ L 2 ^) dependent 
on the nature of the disorder. The Larkin length corre- 
sponds to the length for which the displacements are of 
the order of the size of the interface or the correlation 
length of the random potential [53j B(L C ) = max(£, rt). 
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The Larkin length is thus the smallest length at which 
the wall can be weakly pinned, and above which it can 
adjust elastically to optimize its local configuration [H3 |. 
Above R r one can thus write 



B(L > L c ) 



mnxli,. r r r [ — 



2C 



(8) 



The roughness exponent £ is a function of the type of 
disorder present in the film, and the dimensionality of 
the manifold. For a line [d = 1) and purely thermal 
fluctuations in the absence of disorder, £ = 0.5. In a 
random bond scenario, an exact value of Q = 2/3 has 
been calculated for a line pi Is! IbH] and C ~ 3/5 
l38l | is expected for the two-dimensional manifolds giving 
values of \i = 0.25 and fi ~ 0.5 - 0.6, respectively, for 
the dynamical exponent /i in these scenarios We 
note that random bond disorder exponents have been 
confirmed by measurements of domain wall creep and 
roughness in an ultrathin magnetic film ^3] • In a random 
field scenario, in which defects locally asymmetrize the 
ferroelectric double well potential, ( = [3!j, giving 
fi = 1 for all dimensionalities of the manifold between 1 
and 4. 

L c is also the length scale at which pinning appears in 
the system in the presence of a driving force. Using 
© for u ~ £ and L = L c one can directly obtain|53 the 
critical field E r 



p. 



(9) 



For driving forces above the critical force / c , the interface 
is unpinned even at zero temperatures since the force is 
large enough to overcome the pinning barriers. For / <C 
f c , however, the force is not large enough to overcome the 
barriers and the motion proceeds by thermal activation. 
This is the creep regime, leading to a small and non-linear 
response. In our case, a rough estimate of the values 
of E c may thus be obtained by extrapolating the linear 
behavior of the velocity, which occurs at high field values. 
Although we were unable to extend our measurement 
significantly into this region, we can nonetheless at least 
place a lower bound on the value of E c of 180 MV/m, as 
indicated on Fig. 1111 for one of our thinner films, where 
higher values of the field could be implemented. Taking 
£ to be of the order of a unit cell, we can use the field 
data to extract an approximate value of L c ~ 0.2 nm, 
below the limit of resolution of our measurement. 

In the creep regime, we can rewrite (0. For simplicity 
we write formulas for the isotropic case. Using the scaling 
u ~ ^(L/Lc)^ one obtains 



E(u, L) = U c 



d-2+2C 



2P s EL d f 



d+Q 



(10) 



where U c — a w ^ 2 Lf 2 . Minimizing the energy with re- 
spect to the external field E, we obtain the size of the 




2.0 3.0 
electric field (10 8 V/m) 

FIG. 11: Domain wall velocity as a function of the applied 
electric field in a 37 nm film. Extrapolating the linear behav- 
ior at high fields allows the critical field E c to be estimated 
as 100 MV/m. 



minimal nucleus capable of moving to a lower energy 
state as 



J creep 



/L e = (/ c //) 1/(2 - C) 



(11) 



with / = 2P S E. The minimal barrier height to be passed 
by thermal activation thus corresponds to the length L* , 
leading to a velocity of the form 



v^exp(-(3U c (fc/f)^^) 



(12) 



if one assumes an Arrhenius law in passing the barri- 
ers. The very slow (creep) response is due to the fact 
that for a small force the system would have to rear- 
range large portions of the interface to be able to find a 
new metastable state of low enough energy. The barriers 
a domain wall must pass to make such a rearrangement 
therefore diverge as the force goes to zero. This diver- 
gence could inherently explain the very high stability of 
domain structures we have observed in our films, in the 
limit of zero electric fields, even at relatively high tem- 
peratures. 

The expression <|llfl gives the critical nucleus size 
L creep as a function of the applied field E and L c . We 
note that this expression is independent of the dimen- 
sionality of the film, and that the applied and critical 
fields are present as a ratio, thus removing the uncer- 
tainty associated with the correction of the field in the 
AFM tip-ferroelectric thin film configuration. As for the 
case of the periodic potential, these expressions are valid 
if the size of the nucleus is smaller than the thickness of 
the sample. Otherwise one of the dimensions of the nu- 
cleus should be replaced by the thickness, transforming 
a two-dimensional interface into a one-dimensional line. 
A crucial difference between the periodic and the disor- 
dered cases is that creep due to disorder can still exist 
in the one-dimensional situation, contrary to the peri- 
odic case. Note that the question of whether the films 
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should be considered as one- or two-dimensional depends 
on which mechanism controls the nucleus. A film could 
thus be in the one-dimensional limit for the periodic po- 
tential, thereby invalidating the periodic potential as a 
possible origin for the creep process, and still be in the 
two-dimensional limit for the disorder provided that the 
size of the nucleus due to disorder remains smaller than 
the thickness of the film. Although creep is still present in 
the one-dimensional disordered case, the value of the ex- 
ponent /j, depends on the dimension. Using the values for 
E c and L c we had obtained, we can estimate the size of 
the critical nucleus for the creep process and compare it 
with that found for the Miller- Weinreich formulation. In 
our system, the applied field is a function of the distance 
r away from the tip center. Using the largest possible 
(random field) value of Q we find L creep to vary between 
0.2 and 1 nm in the thinnest films (29.0 - 51.0 nm), and 
0.2 and 2.5 nm in the thickest films (95.0 - 130.0 nm). 
Note that in all cases L creep ~ O.Olr, where r is the radius 
of the domain, so the approximation of a linear domain 
wall is quite reasonable at that scale +[5^|. 

Interestingly, the anomalously large size of the nucleus 
given by the Miller- Weinreich model has previously been 
remarked on by Landauer [2cj- Moreover, studies of 
the piezoelectric effect and dielectric permittivity in PZT 
films 0, 01 have shown non-linear features which can- 
not be described by a simple phenomenological model, 
but which could be described by the pinning of domain 
walls at randomly distributed pinning centers. |59| The 
scenario of weak pinning of elastic interfaces in disor- 
dered media was also invoked as a possible explanation 
for the dielectric dispersion observed in the ferroelec- 
tric RbH2P04 |43j. The present study, investigating the 
static and dynamic behavior of individual ferroelectric 
domains with nanoscale resolution therefore provides an 
important clarification of this issue from the microscopic 
point of view. 



VI. DOMAIN WALL DYNAMICS IN THE 
PRESENCE OF ARTIFICIAL DEFECTS 

Another way to approach the question of ferroelec- 
tric domain wall creep, and its microscopic origins, is 
to investigate the effects of changing the disorder in the 
film. We studied two different types of defects: colum- 
nar tracks of amorphous material introduced by heavy 
ion irradiation (carried out by C. Simon and A. Ruyter 
at GANIL), and planar a-axis oriented inclusions intro- 
duced in thicker films during growth. For the irradiated 
films, samples 51.0 to 116.0 nm thick were half-shielded 
by metallic plates before irradiation, in order to directly 
compare domain dynamics in non-irradiated and irradi- 
ated regions of the same sample. As shown in the AFM 
topography scan in Fig. ll2f ah the tracks are 6-13 nm in 
diameter at the point of entry, and appear as small white 
"bumps" , since the amorphous material is less dense than 
the crystalline form surrounding the tracks. As shown 
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FIG. 12: Topography (a) and PFM (b) measurements of a 
95.4 nm sample with columnar defects created by heavy ion 
(Pb) irradiation. On the film surface, the defects appear as 
white "bumps" in the topography, corresponding to the dark 
contrast regions in the PFM image, which we were unable to 
switch to the uniform polarization of the background. 

in the corresponding piezoresponse image (Fig. I12f b)h 
these are exactly the regions which remain unaffected by 
attempts to uniformly polarize the area by scanning with 
a constant negative voltage, appearing as dots with dark 
phase contrast, as opposed to the lighter phase contrast 
of the rest of the image. Thus, irradiation defects can be 
distinguished from merely topographic features like the 
large white "bump" in the middle of the topographical 
scan, which switches under the application of an electric 
field and does not appear as a corresponding dark area in 
the piezoresponse image. The density of defects observed 
in this image (4.7 x 10 10 /cm 2 ) agrees reasonably with 
the nominal irradiation density 4 x 10 10 /cm 2 used, con- 
sidering the small area (300 x 300 nm 2 ) over which the 
Poisson distribution is observed. To ensure continuous 
tracks throughout the sample, with the ions lodging deep 
in the substrate, measurements from other experiments 
on similar oxides carried out at GANIL were used to de- 
termine irradiation energies |44|. To produce the planar 
a-axis inclusion, we grew relatively thick films, slightly 
varying the temperature from its optimal range, which 
generally results in the presence of a-axis, as shown in the 
AFM topography (Fig. I13f a)h Once again, the a-axis 
oriented regions do not respond to switching attempts, 
as shown in the piezoresponse image in Fig. I13f . since 
the applied field is perpendicular to the polarization axis. 
These defects are approximately 25 nm wide, extending 
for hundreds of nanometers in the crystallographic a and 
b axis directions, and inclined at 45° to the c-axis of the 
film |45|. In all films, 25 domain arrays were written, and 
the domain dynamics extracted as described previously 
|27|. In the a-axis films, domains were written adjacent 
to, but not right on top of the a-axis inclusions. 

As shown in Fig. ll4f a') for a 95.4 nm film, we find that 
domain sizes are comparable in both the irradiated and 
normal regions of the sample. However, for longer writ- 
ing times, domains in the irradiated part of the sample 
appear somewhat larger that those written in the nor- 
mal regions of the same sample. When comparing pure 
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FIG. 13: Topography (a) and PFM (b) measurements of a 
108.9 nm sample with a-axis inclusions. These appear as 
criss-crossing lines on the surface, parallel to the a and b 
crystallographic axes of the film. Since the direction of the 
polarization is in plane, these regions do not switch when a 
perpendicular electric field E is applied across the film, once 
again appearing as dark contrast lines against a uniformly 
polarized background. 



c-axis film 95.4 and 130.0 nm thick to a film of similar 
thickness (108.9 nm) with a-axis inclusions (Tig. EP 3 )) 
this increase in domain size for longer writing times in 
the presence of defects is also observed. The effect of the 
irradiation tracks and a-axis inclusions becomes more ev- 
ident when the Arrhenius plot of velocity vs the inverse 
field is compared for the same films fFig. H5f a1 and (b)). 
We find a decrease in the values of the dynamical expo- 
nent in the irradiated regions, as compared to the normal 
regions, from a range of 0.62 - 0.69 to one of 0.38 - 0.5. 
This decrease is even more marked for films containing a- 
axis inclusions, with \i values of 0.31 and 0.19 [2f^ . That 
this is not a thickness effect can be seen by comparing 
the a-axis-containing films to pure c-axis films of similar 
thickness, in which fi values of 0.69 and 0.78 (for 95.4 
and 130.0 nm thick films, respectively) are observed. It 
is interesting to note that both these types of defects 
are large compared to the thickness of the wall itself, 
occur at a relatively low density, and are either colum- 
nar or planar, as opposed to point-like. The mechanism 
by which they affect the dynamics of the domain wall is 
therefore unlikely to be a direct weak collective pinning 
effect, and thus the marked decrease observed in the dy- 
namical exponent is somewhat surprising. Possibly, the 
relaxation of strain in the film caused by these macro- 
scopic "defects" changes the density or pinning force of 
the inherent, weak point defects present in the film, thus 
affecting the disorder potential experienced by the do- 
main wall. Since a large amount of energy is dissipated 
within the substrate by the heavy ion during irradiation, 
with possible deformation as strain induced in the sub- 
strate [46J, |47j , this could result in additional strain ef- 
fects on the ferroelectric film. Alternatively, the presence 
of defects which penetrate through the entire thickness 
of the film might act to make the wall more rigid - and 
therefore more one-dimensional. A more detailed inves- 
tigation of domain walls in the presence of these defects, 



1 

a 

• i 

a 

o 

T3 



150 
100 
50 







V normal 
T irradiated 



954 V ft 



V* 1 



(a) 



10" 5 10" 3 10" 1 10 1 
writing time (s) 



200 



t/3 

3 

'S 

S-H 

a 

£ 

o 



100 







v normal, 954 A ▼ 
A normal, 1300 A tT t 
t a-axis, 1089 A T 



(b) 



ia 6 ia 4 i0" 2 i io 2 

writing time (s) 



FIG. 14: Domain radius as a function of the writing time 
for the 95.4 nm half-irradiated film (a), and for the 108.9 nm 
film with a-axis inclusion (b). In both cases, the domains 
written with longer writing times appear larger than in the 
non-irradiated or pure c-axis films included as a comparison. 



perhaps looking at their specific individual interaction, 
is needed to ascertain the microscopic nature of this be- 
havior. 



VII. CONCLUSION 

Using the unprecedented control and precision pro- 
vided by AFM, we were able to study the growth of 
individual nanoscale ferroelectric domains in epitaxial 
thin films, investigating both the fundamental physics 
of domain wall motion, and its consequences for possible 
AFM-ferroelectric applications. Our studies demonstrate 
that domain wall motion in ferroelectric thin films is a 
creep process in which v oc exp(— [3U C {E C / E)* 1 ), with a 
dynamical exponent ^ between ~ 0.7 and 1.0. This pro- 
cess controls the lateral growth of domains in low electric 
fields applied by an AFM tip. A detailed analysis of the 
possible microscopic origins of the observed domain wall 
creep suggests that it is the result of competition between 
elastic behavior and pinning in a disorder potential. The 
reduced dimensionality of our thin films compared to the 
size of the critical nucleus precludes pinning in the com- 
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FIG. 15: Domain wall speed as a function of the inverse ap- 
plied field for the films shown in Fig. 1141 In the presence of 
defects, the values of the creep exponent fj, decrease from 0.69 
to 0.38 and from 0.69 and 0.78 to 0.31 in the half-irradiated 
film, and the pure c-axis films vs the film with a-axis inclu- 
sions, respectively. 



mensurate potential of the crystal itself as the mecha- 
nism for the non-linear field dependence of the velocity. 
We have also observed that the presence of artificially in- 
troduced defects, such as irradiation tracks or a-axis in- 
clusions, strongly decreases the dynamical exponent for 
domain wall creep, from ~ 0.6 - 0.7 to 0.4 - 0.5 and 
~ 0.31. All the domains show high stability (up to 4 
months for the longest duration experiments), inherently 
explained by the physics of a system in which elasticity 
and pinning by a disorder potential compete, leading to 
glassy behavior in the presence of low electric fields. This 
high stability is important for applications such as infor- 
mation storage or surface acoustic wave devices based on 
ferroelectric domain structures. Our studies have also 
identified the key parameters controlling domain size as 
the writing time, the applied voltage, and the confine- 
ment of the electric field, with encouraging results for 
novel tip geometries using carbon nanotube bundles. 
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[57] As before the length here is the length perpendicular to 

the polarization direction. 
[58] In the thinnest films, longer writing times resulted in un- 



stable and destructive interactions between the tip and 
the sample. Therefore the larger domains obtained for 
long writing times were only measured in the thicker 
films. 

[59] We note that the studies referred to were carried out in 
thick (over 1 /im) ceramic and sol-gel films, where the 
presence of multiple grain boundaries and differently ori- 
ented domain walls provides a much more complex dis- 
order landscape compared to the epitaxially grown films 
used for this study. 



